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CALCULATION  OF  THE  BIOT-SAVART  LAW  MAGNETIC  FIELD 
USING  CURRENT  DISTRIBUTIONS  OBTAINED 
FROM  A  FINITE  ELEMENT  ANALYSIS 


INTRODUCTION 

Three-dimensional  static  magnetic  field  solutions  using  finite 
elements  have  been  formulated  usinq  both  scalar  and  vector  potential 
functions.  Guancial  and  DasGupta,  Frye  and  Kasper,2  and  Zienkiewicz3 
have  addressed  the  vector  potential  formulation. 

The  use  of  scalar  potential  functions  provides  great  computational 
advantages  over  vector  potential  functions,  since  the  number  of  unknowns 
at  a  grid  point  is  reduced  from  three  to  one.  Zienkiewicz  et  al.4 
showed  that  a  finite  element  scalar  formulation  was  possible  using  the 
Biot-Savart  law  to  account  for  current  terms.  Armstrong  et  al.5  indi¬ 
cated  how  the  Zienkiewicz  formulation  could  be  improved  to  eliminate 
ill-formed  matrices  by  a  breakup  of  the  problem  region  into  domains 
using  different  scalar  potentials 

Although  the  use  of  the  scalar  potential  function  gives  great  economy 
in  the  solution  of  the  finite  element  equations,  there  are  problems  where 
the  evaluation  of  the  Biot-Savart  magnetic  field  requires  much  computation. 
Typically  such  problems  arise  where  an  object  produces  electric  currents 
internally  in  its  own  structure  and  externally  in  the  infinite  conducting 
medium  in  which  it  is  embedded.  Figure  1  shows  a  diagram  of  such  a 
problem.  Wikswo6  shows  how  the  evaluation  of  the  Biot-Savart  law  magnetic 
field  for  such  a  problem  can  be  reduced  to  an  integration  over  boundaries 
where  there  is  a  change  in  media  conductivity  and  over  regions  where 
current  sources  are  present.  Often,  when  the  magnetic  field  around  such 
an  object  must  be  evaluated,  the  currents  due  to  the  electric  field  are 
calculated  using  finite  element  techniques.  This  is  particularly  true  if 
the  object  is  of  complex  shape  or  has  unusual  electric  field  boundary 
conditions  with  the  infinite  medium.  Typically,  the  finite  element 
electric  field  model  of  an  object  will  be  made  up  of  rod  and  plate 
elements,  although  volume  elements  may  sometimes  also  be  necessary.  The 
infinite  medium  is  modeled  as  a  set  of  volume  elements  with  an  appropriate 
boundary  condition  to  terminate  the  finite  element  mesh.  The  purpose  of 
this  report  is  to  show  the  derivation  of  equations  from  Wikswo 's  formula¬ 
tion  for  evaluating  the  Biot-Savart  law  magnetic  field  by  using  currents 
obtained  from  a  finite  element  electric  field  analysis. 
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THEORETICAL  BACKGROUND 


In  the  scalar  potential  finite  element  formulations  for  static  mag¬ 
netic  fields,4  the  total  field  is  assumed  to  be  equal  to  the  sum  of  two 
components. 


a  =  i  +  a  . 

c  m 

Component  tf  is  the  Biot-Savart  law 
component  Hm  is  the  gradient  of  a  scalar 


(1) 

magnetic  field  strength,  and 
potential  function. 


ac(*> 


=  [ 

4 t r  J.., 


3  x  (3-fr) 

\t-t'  | 3 


dV 


(2) 


V ?*  •  (3j 

where  J  is  some  known  current  density  distribution;  ft  is  the  position 
vector  to  the  point  where  Hc  is  to  be  evaluated,  a*  is  the  position 
vector  of  dV',  the  element  volume  of  the  conductor  with  current  J;  and 
<f>  is  the  unknown  scalar  potential  function. 


Use  of  this  assumed  solution  for  the  magnetic  field  strength 
automatically  satisfies  the  field  relation  between  magnetic  field 
strength  and  current  density: 

v  X  ft  =  ft  X  ac  =  3  (4) 
since  ft  x  ftm  =  ft  x  V4>  =  0  .  (5) 


The  remaining  field  relation  requiring  that  the  divergence  of  the 
flux  density  be  zero  establishes  the  required  relation  for  determining 
the  scalar  potential  function; 

ft*ft  =  0  (6) 

ft*v(ft<l>)  =  -ft*yffc  .  (7) 


The  finite  element  formulation  given  by  Zienkiewicz  et  al.4  for 
the  solution  of  the  field  equation  results  in  a  matrix  equation  that 
relates  a  set  of  discrete  scalar  potential  unknowns  {<(>}  to  a  set  of 
loadings  {f}  obtained  from  the  Biot-Savart  law  magnetic  field  strengths 
via  a  matrix  [Kj  calculated  with  finite  element  interpolation  functions 


[K]  {<}>}  =  {f}  (8) 

Kij  =  /v  ^Ni>*^Nj)  dV  (9) 


* 


*3 
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(10) 


It  is  not  necessary  to  discuss  details  of  the  finite  element 
technique  since  all  such  details  are  standard  and  are  discussed  fully 
elsewhere. 7 

The  solution  of  the  field  problem  proceeds  by,  first,  constructing 
an  appropriate  finite  element  model;  second,  calculating  the  Biot-Savart 
magnetic  field  strengths  from  the  known  current  distributions;  third, 
computing  the  loadings  f.  from  equation  (10);  fourth,  solving  the  matrix 
equation  (8)  for  {$};  and  finally,  fifth,  adding  the  potential  gradients 
obtained  from  the  finite  element  solution  to  the  Biot-Savart  field  strength 
to  obtain  the  total  field  result.  The  finite  element  solution  acts  as  a 
correction  to  the  Biot-Savart  law  field  strength  to  account  for  magnetiza¬ 
tion  existing  in  materials  with  permeability  other  than  that  of  free  space. 

For  problems  with  infinite  conductor  domains,  the  evaluation  of  the 
Biot-Savart  law  magnetic  field  value  from  equation  (2)  requires  a  great 
amount  of  computation  because  the  diffuse  nature  of  the  current  densities 
results  in  an  extensive  required  region  of  integration.  Wikswo6  has  shown, 
in  the  case  where  current  densities  are  known  to  approach  zero  in  regions 
of  the  infinite  conducting  medium  removed  from  the  domain  of  interest,  that 
equation  (2)  can  be  replaced  by  the  following  relation: 


V'  x  3(ft') 

I*  -  *  I 


dV’. 


(ID 


The  advantage  to  this  relation  is  that,  since  steady  currents  can  be 
related  to  the  gradient  of  a  scalar  potential  function,  the  curl  of  the 
current  density  is  zero  at  all  points  in  the  volume  of  the  problem  except 
at  surfaces  where  there  is  a  change  in  conductivity  or  a  current  source. 
Thus,  the  integration  of  equation  (11)  reduces  from  a  volume  integration 
to  a  surface  integration.  For  practical  purposes  we  can  break  this  sur¬ 
face  integration  into  that  over  the  surface  of  rods,  plates,  and  volume 
conductivity  interfaces.  Most  metal  structures  can  be  characterized, 
for  the  purpose  of  defining  current  distributions,  as  a  collection  of 
rods  that  carry  currents  in  a  one-dimensional  manner  and  of  plates  with 
two-dimensional  current  sheets.  The  currents  in  the  conductive  medium 
in  which  the  metal  body  is  embedded  can  be  accounted  for  by  integrating 
over  the  interface  between  the  surface  of  the  metal  body  and  the  conductive 
medium,  and  along  the  interfaces  of  differing  conductivity  in  the  media. 

In  the  following  sections,  separate  equations  will  be  developed  from 
equation  (11)  for  the  rod  and  plate  surfaces  and  the  volume  interface 
surfaces.  The  equations  are  developed  based  on  the  assumption  that  the 
current  density  outside  the  region  of  the  rod,  plate,  or  volume  is  zero. 


£ 
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It  will  later  be  shown  that  the  summation  of  separate  integrations  appro¬ 
priately  accounts  for  the  situations  where  this  assumption  is  not  the 
case. 


INTEGRATION  FOR  RODS 


Figure  2  shows  the  rod  geometry  and  the  coordinate  systems  used  to 
describe  the  rod  integrations.  The  rod  axis  is  assumed  to  lie  along  a  line 
between  two  points  in  space  called  grid  points  1  and  2.  The  two  grid  points 
have  their  position  in  space  defined  by  global  coordinates  x,  y,  and  z.  The 
rod  axial  coordinate  ZR  is  taken  as  having  its  origin  at  grid  point  1  and 
as  having  a  positive  sign  in  the  direction  of  grid  point  2.  The  axial 
coordinate  direction  vector  is  given  by  zR. 


z  =  *  1*3-.-  *-0  +  p  iy.2_:  lA  +  z  z.d  . 

|t*2  -  $i|  |$2  ~  3l|  |$2  ~  %i\  (12) 

The  rod  coordinate  direction  vector  xr  is  taken  as  the  cross  product 
of  the  global  y  direction  and  the  zr  direction*. 

/V  A 

X  y  x  ZR 

^  X  ZR'  (13) 

In  the  case  where  y  and  zr  are  parallel,  XR^is  taken  as  being  the  same 
as  xA  The  rod  coordinate^ rection  vector  yR  is  the  cross  product  of 
the  zr  direction  and  the  xr  direction; 

yR  =  ZR  X  XR  *  (14) 

The  xr  and  9r  direction  vectors  can  be  expressed  in  terms  of  their 
respective  components  as 


R 

>R 


Rx  +  y  XRy  +  2  XRz 

(15) 

Rx  +  *  YRy  +  *  YRz  • 

(16) 

The  rod  is  assumed  to  have  a  circular  cross  section  and  to  have  a 
radius  (a)  that  is  small  with  respect  to  the  distance  to  the  point  where 
the  Biot-Savart  law  magnetic  field  value  is  to  be  calculated.  In  other 
words,  the  calculations  are  not  to  be  made  in  the  nearfield  of  the  rod, 
at  least  from  the  standpoint  of  rod  cross-sectional  geometry. 
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a  «  |ff  -  ff'  |  (17) 

*  Using  cylindrical  coordinates,  we  can  define  the  radial  rod  direction 
r  and  tangential  direction  t  in  terms  of  the  xR  and  yR  directions  as 


r  = 

cos0  xR  +  sin0  yR 

(18) 

/N 

t  = 

-sin©  xR  +  cos©  yR  . 

(19) 

The  vector  ff  -  ff'  can  be  rewritten  as 

ff  - 

ff'  =  [ff  -  (ffc  +  r)]  , 

(20) 

where  ff„  : 

c 

=  X  (Xi  +  zR*x)  +  y  (yi  +  zR*y)  +  z  (z1  +  zR.z) 

(21) 

ZR  : 

=  ZR  ZR 

(22) 

-V 

r  = 

A 

r  r  . 

(23) 

zR  and  r  are  the  amplitudes  of  the  axial  and  radial  coordinates, 
ffc  is  the  vector  to  a  point  on  the  rod  axis.  ^  can  also  be  written  in 
terms  of  its  global  coordinates  as 

ffc  =  V  +  V  +  v  •  (24) 


Substituting  equations  (15),  (16),  (18),  (19),  and  (24)  into  equation 
(20)  and  evaluating  th|  absolute  value  of  the  vector  expression  gives  the 
following  result  for  |ff  -  ff' | : 


-  ff'|  =  |(x  -  xc  -  r 

+  (y  -  yc  -  r 

+  (Z  -  zc  -  r 

Since  r  will  always  be  small 
be  accurately  approximated  as 


LcoseX^  +  sineYRx])2 
[coseX^y  +  sineY^])2 
[cos0XRz  +  sineYRz])2}%  . 
in  relation  to  |ff  -  ffc|»  |ff 


(25) 
ff'  |  can 


Iff  -  ff'| 


-  r(x-xc)  [cos0XRx  +  sinoY^] 


-  r{y-yc)  [cos0XRy  +  sin0YRy] 

-  r(z-zc)  [cos0XRz  +  sin0YRz;A  •  (26) 
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The  current  In  the  rod  Is  assumed  to  be  flowing  in  the  axial  direction 
of  the  rod;  thus,  Jrx  and  Jrw  components  of  the  J  current  density  are  zero. 
The  curl  of  J  in  terms  of  roa  coordinates  is  given  by 


A  /V 

r  t 


/s 


V 


'  x  J  = 


_L  JL  JL 

3r  at  azR 


0  0JRz 


[%]  -  [-  ^] 


+  (0).  (27) 


3J 


Rz 


3t 


If  the  rod  current  density  is 
is  always  zero, even  along  the 


Thus,  the  curl  of  J  simplifies  to 
v  x  J  =  -t 


Rz 

3r 


assumed  constant  over  the  rod,chen 
sides  of  the  rod  parallel  to  the  axis. 


(28) 


9^rz 

The  derivative  ■■■■■■  is  zero  at  all  points  on  the  rod  except  along 

the  rod  longitudinal  surface, where  the  value  of  Jrz  is  assumed  to  drop 
to  a  zero  value. 


In  evaluating  the  integral  of  equation  (11)  it  is  necessary  to 
evaluate  only  the  limit  over  a  film  of  volume  at  the  rod  longitudinal  sur¬ 
face  as  the  film  thickness  goes  to  zero.  Figure  3  illustrates  such  a 
film  of  volume.  The  differential  volume  dV'  can  be  expressed  as 


dV'  =  a  6r  de  dz, 


(29) 


9^rz 

The  partial  derivative  r-  can  also  be  expressed  as  the  following 
limit:  °r 


8JRz  _  lim  JRz(a  +  T~)  -  JRz(a  ‘  T") 


3r  6r->o  6r 

The  combination  of  the  terms  of  equations  (29)  and  (30)  produces 
the  following  expression; 

5r>  ,  6r 


(30) 


9J 


Rz 


3r 


dV 


=  lim  JRzU  +  ~  JRz^a  ~  X! 


r-*-a 


6r-K) 


6r 


a  ar  de  dzR  .  (31) 


JRz  (a  +  4r)is  interpreted  to  be  the  current  density  just  outside 
the  rod  volume  and  is  zero.  The  current  density  JRz  (a  -  4£)  is  the 


8 


ROD  SURFACE 


Figure  3.  Evaluating  the  Limit  of  a  Film  of  Volume 
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current  density  just  inside  the  rod  surface  and  is  equal  to  J-  .  Thus, 
expression  (31)  becomes, after  the  evaluation  of  the  limit. 


■  -JR2ad6d2R  • 


The  volume  integration  for  the  rod  given  by  equation  (11)  then 
reduces  to  the  following  integration  over  the  rod's  longitudinal  surface: 


L  2tu 


nc  <*) 


.±  f  fll 

4tt  J  I R- 


adedz. 


The  integration  limit  from  0  to  2tt  can  be  changed  to  a  limit  that 
varies  from  0  to  tt  by  adding  up  the  contributions  of  elemental  surface 
areas  at  0  and  e  +  tt  radians  simultaneously: 


The  first  term  in  the  integration  accounts  for  the  elemental  area 
at  e,and  the  second  term  accounts  for  the  elemental  term  at  0  +  ir.  Here 
we  have  accounted  with  a  minus  sign  for  the  fact  that  the  tangent  vector 
at  9  +  1  is  in  the  opposite  direction  from  the  tangent  vector  at  0.The 
terms  |R-R' L  and  | R-R 1 |0+7r  can  be  calculated  from  equation (26)  by  setting 
r  =  a;  0 


If  equations  (35)  and  (36)  are  substituted  into  equation  (34)  and 
the  expression  inside  the  integration  is  placed  under  a  common  denominator, 
we  obtain  the  following: 


10 
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Performance  of  the  integration  with  respect  to  6,  which  involves  the 
algebraic  expansion  of  equation  (37),  results  in 


ac  <«> 


■WH1  x  X  >R 
iHi3 


(38) 


If  Id  is  the  total  current  flowing  in  the  rod,  then  Jrz  is  equal  to 
Ir  divided  by  the  rod  cross-sectional  area: 


JRz  "  ia2-  • 

AAA 

Also  zR  =  xR  x  yR 


(39) 

(40) 


Substitution  of  equations  (39)  and  (40)  into  equation  (38)  results  in 
the  final  form  for  the  integration  of  the  rod  current  contribution  to  the 
Biot-Savart  law  magnetic  field  strengths 


ifc  (S) 


IRZR  X 

iHi3 


(41) 


The  integration 
then  reduced  for  the 


with  respect  to  a  volume  given  by  equation  (11) 
rod  to  an  integration  along  the  rod  axis. 


is 


INTEGRATION  FOR  PLATES 


Figure  4  shows  the  coordinate  systems  used  in  the  development  of 
expressions  for  integrating  over  the  region  of  plates.  The  plate  shown  in 
the  figure  is  quadrilateral,  but  the  coordinate  system  definition  is  essen¬ 
tially  the  same  for  triangular  plates.  The  plate  is  defined  by  four  grid 
points,  at  the  vertices  of  the  plate  edges,  that  have  their  locations 
specified  in  terms  of  global  x,  y,  and  z  coordinate  directions.  These 
points  are  known  as  grid  points  1,  2,  3,  and  4.  The  normal  to  the  plate 
surface  can  be  calculated  by  taking  the  cross  product  of  a  vector  directed 
between  the  first  and  second  grid  points  and  a  vector  directed  between  the 
first  and  fourth  grid  points  and  normalizing  the  result: 


a  x  $ 
|cx  x  $| 


(42) 
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Figure  4.  Coordinate  Systems  for  Integrating  Over  the  Region  of  Plates 


where  a  =  (x2  -  Xj)  x  +  (y2  -  yx)  y  +  U2  -  Zj)  z  (43) 

I  ■  (x4  -  Xj)  x  +  (y4  -  yx)  y  +  (z4  -  Zj)  z  .  (44) 

xl»  yi.  zl»  x2»  n*  z2»  an<*  x4»  ^4»  z4»  are  9lobal  coordinates  of 
the  first,  second,  and  fourth  grid  points, respectively. 

A 

The  first  tangential  plate  coordinate  direction, known  as  ti»is  in 
the  same  direction  as  a  ; 

‘i  ■  iSr  •  (45) 

The  second  tangential  coordinate  direction  is  obtained  by  taking  the 
cross  product  of  n  and  * 

t2  =  n  x  tj  .  (46) 

The  advantage  of  using  these  coordinate  directions  is  that  finite 
element  programs  that  analyze  current  fields  often  use  these  coordinates 
when  listing  components  of  a  current  vector. 

The  curl  of  the  current  density  for  a  current  flowing  in  the  plate 


^  ov2J  ^  l  °  J  1  ^  J 

Note  that  there  is  no  component  of  current  density  normal  to  the 
plate  as  all  currents  are  assumed  to  be  flowing  parallel  to  the  plane  of 
the  plate. 

The  current  in  the  plate  can  be  expressed  as  being  equal  to  the 
plate  material  conductivity  times  the  gradient  of  a  potential  function; 


Jt2  ■ « %  •  (49> 

Substituting  equations  (48)  and  (49)  into  the  first  term  of  equation 
(47),  we  have 


I  =  _ ?.2±.  =  o  . 

ni.  m.  A 1  w 


3tj8t2 


8t13t2 


The  component  of  the  curl  term  normal  to  the  plate  drops  out.  The 
remaining  terms  in  the  curl  relation  are  normal  derivatives  of  current 
terms  tangent  to  the  plane  of  the  plate.  These  terms  are  zero  everywhere 
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except  at  the  surface  since  the  current  is  assumed  not  to  vary  across 
the  plate  thickness.  At  the  surface,  the  current  vector  amplitude  changes 
abruptly  from  some  value  to  zero.  Figure  5  shows  a  diagram  of  the  ele¬ 
mental  volumes  at  the  upper  and  lower  surfaces  of  the  plate.  The  expres¬ 
sion  for  the  elemental  volume  may  be  written  as 

dV  =  STdt1dt2.  (51) 


30  f  3J. 

l 

The  derivatives  and  can  also  be  written  as  limits 
an  an 

8Jt  k  (r  +  T-}-Jt 

_h  =  lim  tl  1  2  t1  2  2 

3n  ~  6T>o  6T 


)Jt  Jt  (I  +  '  Jt  (!  ‘  T> 

Z2  =  lim  z2  *  d  h  2  2 

3n  6T-*o  L  5T 


On  the  upper  plate  surface  the  current  densities  above  the  plate 
are  zero,  and  on  the  lower  surface  the  current  densities  below  the  plate 
are  zero: 

UPPER  SURFACE  0t  (£  +  =  J  (|  +  ^)  =  0  (54) 

LOWER  SURFACE  J  (^  -  ^)  =  Jt  (j  -  ^)  =  0  .  (55) 


On  the  upper  surface  the  current  densities  below  the  surface  are 
the  plate  current  densities,  and  on  the  lower  surface  the  current  den¬ 
sities  above  the  surface  are  the  plate  densities; 


UPPER  SURFACE 


Jt1  (jr  "  Tf)  =  Jtj 
Jt2  ^7  '  =  Jt2 


LOWER  SURFACE 


\  (?  +  &  -  Jt, 
Jt2  +  =  Jt2 
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Substituting  equation  (47)  and  equations  (51)  through  (59)  into 
equation  (11)  gives 


where  |1$  -  ft'|  is  the  amplitude  of  the  vector  to  the  upper  surface 

and  -  ^'|-|s  1S  the  amplitude  of  the  vector  to  the  lower  surface  . 

The  amplitudes  |$  -  lus  and  -  ft' (is  can  be  written  in  terms  of 
vectors  to  the  plate  midsurface  and  the  normal  of  the  plate; 

l*-*'lu$-  ls-Vr"l  (6l) 

I*  -  *'lis "  Is  *  •  (62) 

where  1L  is  the  vector  to  the  midsurface. 


Under  the  assumption  that  the  plate  thickness  is  small  in 
to  the  distance  to  the  point  where  the  field  strength  is  to  be 
equations  (61)  and  (62)  can  be  approximated  by 


Is  •  S'lusK1j7jy  jj*  -  \,!2  *^"(S  '  v] 


relation 
calculated , 


(63) 


I*  -  *  I 


Is 


is  \ 


(64) 


Substituting  equations  (63)  and  (64)  into  equation  (60)  and  placing 
the  two  terms  in  the  integral  under  a  common  denominator  results  in 


If  the  sheet  current  ts  is  defined  as  the  current  density  times  the 
thickness  of  the  plate,  then  equation  (65)  can  be  rewritten  in  the  follow¬ 
ing  way; 


Thus, the  Integration  of  the  plate  Is  reduced  to  an  Integral  over 
the  midsurface  of  the  plate. 
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VOLUME  INTERSECTION  SURFACE 

Figure  6  shows  the  coordinate  systems  used  in  the  development  of 
expressions  for  integrating  over  surface  regions  of  volumes.  The  region 
shown  in  the  figure  is  quadrilateral,  but  the  coordinate  system  definition 
is  the  same  for  triangular  regions.  The  surface  region  is  defined  by  four 
grid  points,  at  the  vertices  of  the  region's  edges,  that  have  their  loca¬ 
tions  specified  in  global  x,  y,  and  z  coordinate  directions.  The  grid  points 
are  numbered  counterclockwise  when  viewed  from  outside  the  conductive  medium. 
The  vector  directions  n,  ti,  and  tp  are  defined  in  the  same  manner  as  with 
the  plate  coordinates,  i.e.,  with  equations  (42),  (45),  and  (46). 


In  the  case  of  a  surface  region  of  a  volume, there  may  be  components 
of  current  density  flowing  in  a  direction  normal  to  the  surface  as  well 
as  tangent  to  it.  Thus. the  curl  of  the  current  density  contains  all  com¬ 
ponents  of  the  current  density: 


As  with  the  plate  element, the  normal  component  of  the  curl  of  3 
always  vanishes  because  the  current  density  can  be  written  as  the  gradient 
of  a  scalar  potential  function  in  the  region  of  the  conductive  medium. 

In  fact, all  components  of  the  curl  vanish  at  all  points  in  the  domain  of 
the  conductive  medium  except  at  points  on  the  surface  region  where  the 
current  density  vanishes  just  outside  the  surface. 

Figure  7  shows  a  diagram  of  the  elemental  volume  of  integration  on 
the  surface  region  of  a  volume.  The  expression  for  the  elemental  volume 
may  be  written  as 

dV  =  Sudtjdtg  .  (68) 


The  derivatives 


-—can  be  written  as  the  following  limits: 


(69) 

(70) 
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Jt  (4^-)  and  Jt  refer  to  current  densities  just  outside  the 

conductive  medium  and  are  both  zero.  Jt  (-  4^-)  and  Jt  (-  4r)  refer  to 


current  densities  just  inside  the  conductive  medium  and  are  equal  to 
J.  and  J.  .respectively. 

h  z2 

Using  equations  (67),  (68),  (69),  and  (70)  in  equation  (11)  and  noting 
that  the  curl  term  in  the  direction  normal  to  the  surface  vanishes,  we 
obtain 


ac«)  - 


A 

h 

lim  | 
Su-K) 

v 

y  atj  +  WJ 

U] 

i  iiJ  _  i*  i 

±[ 

4it  / 

S  . 

r  / 

t  LJl.Jl 

2  |j5u-K)  y  3tj  <5u 


Su 


dt^dt2 


(71) 


Evaluation  of  the  limits  of  equation  (71)  results  in  the  final  expres¬ 
sion  for  the  integration  over  the  surface  region.  In  evaluating  the  limit 
it  should  be  noted  that  the  tangential  derivatives  of  the  normal  component 
of  current  are  finite  since  the  discontinuity  in  the  normal  component  of 
current  density, if  it  exists, is  in  the  direction  normal  to  the  surface, 
not  tangential  to  it: 


Bc<*)  ■  i  / 


x  n 


I*  -  *'l 


ds 


(72) 


It  is  interesting  that  the  components  of  the  current  density  normal 
to  the  surface  region  do  not  enter  into  the  magnetic  field  strength  cal¬ 
culations.  Only  tangential  current  density  components  are  of  importance. 


SURFACES  WITH  NONZERO  CURRENT  DENSITIES  ON  BOTH  SIDES 

In  all  of  the  discussions  thus  far, integration  equations  have  been 
developed  for  surfaces  that  have  current  densities  on  one  side  of  the 
surface  and  not  the  other.  To  consider  the  case  where  a  current  density 
exists  on  both  sides  of  a  surface, it  is  necessary  only  to  add  the  results 
of  two  separate  integrations  where  the  surface  normal  is  defined  separately 
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and  with  a  different  positive  direction  in  each  integration.  Consider 
figure  8, which  shows  a  surface  r  between  two  regions,  A  and  B,  of  dif¬ 
fering  conductivity.  Ja  is  the  current  density  in  region  A,  and  3b  is 
the  current  density  in  region  B. 

In  the  evaluation  of  the  Biot-Savart  law  magnetic  field  strength 

we  have  seen  that  it  is  necessary  to  evaluate  the  integral  /  — — i-  dV  over 

the  surface  of  the  intersection  of  two  regions  of  differing  conductivity. 
To  do  this,  the  limit  is  taken  of  the  product  of  the  derivative  and  the 
infinitesimal  volume  element  as  the  side  of  the  element  approaches  zero*. 


JT  -  JT 

_!a _ |b.  .n 

6ng  6nB 


dS  • 


(73) 


Here  we  have  taken  6n0  as  positive-pointing  away  from  the  B  region. 
Equation  (73)  can  easily  be  rewritten  as 


region  outside 
A  has  zero  cur¬ 
rent  density 


region  outside 
B  has  zero  cur¬ 
rent  density  (74) 
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REGION  OF 
INTEGRATION 


REGION  B 

. 

Jfi  =  aBW 


REGION  A 
JA  *  °P$4> 


SURFACE  OF 
INTERSECTION  BETWEEN 
REGIONS  A  AND  B 


Figure  8.  Surface  Between  Two  Regions  of  Differing  Conductivity 
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The  integration  over  the  surface  of  the  intersection  between  two 
volumes  of  differing  conductivity  is  seen  to  be  equivalent  to  the  sum  of 
the  separate  integrals  over  the  two  regions,  with  the  assumption  that  the 
current  density  in  all  other  adjacent  regions  is  zero.  As  a  result,  the 
Biot-Savart  law  field  strength  may  be  evaluated  by  summing  the  rod,  plate, 
and  volume  intersection  surface  integrations  of  equations  (41),  (66),  and 


This  conclusion  perhaps  could  have  been  reached  in  a  more  direct 
manner  by  simply  noting  that  equation  (11)  is  the  volume  integral  taken  over 
regions  of  nonzero  current  density  and  that  this  integration  can  be  evaluated 
by  summing  the  integration  over  separate  regions. 
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